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ABSTRACT 

Disc spectra from highly luminous black hole binaries are observed to be rather simple, 
despite theoretical predictions to the contrary. We collate the disc-dominated spectra from 
multiple observations of 10 separate sources, and show that these overwhelmingly follow a 
multi-temperature disc blackbody spectrum with luminosity oc T 4 , and that such discs are sta- 
ble. These results are in conflict with standard Shakura-Sunyaev a-discs predictions, including 
proposed modifications such as additional energy loss in a jet, wind or corona. Accretion disc 
spectra can be useful, reliable guides to the disc structure, and should be used to test the next 
generation of accretion disc models, in which the disc viscosity is calculated self-consistently 
from the magnetically generated turbulent dynamo. 
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1 INTRODUCTION 



Accretion onto a compact object is believed to be the source of 
power in X-ray binaries (see e.g. Tanaka & Lewin 1995). At accre- 
tion rates in excess of a few percent of Eddington rate most of the 
black hole binaries dwell in the so-called soft, or high, X-ray spec- 
tral state. Their X-ray spectra are then dominated by a soft thermal 
component, of typical temperature < 1 keV. The standard model 
of accretion (Shakura & Sunyaev 1973, hereafter SS73) explains 
these spectra as emission from an optically thick and geometri- 
cally thin disc. This can form a stable emitting structure down to 
the minimum stable orbit (which depends on the black hole spin), 
so its spectrum contains the imprint of strong special and general 
relativistic effects (e.g. Cunningham 1975; Ebisawa, Mitsuda & 
Hanawa 1991). 

This all makes disc spectra an exciting potential probe of the 
dramatically curved space-time in the vicinity of the black hole. 
However, there are many problems in interpreting them. Firstly, the 
quasi-thermal component which is identified as the disc emission 
is generally accompanied by a power-law tail out to much higher 
energies (e.g. Tanaka & Lewin 1995). This tail is well fit by Comp- 
tonization of seed photons from the disc, requiring that some of 
the energy is not thermalized in the disc but is instead dissipated 
in an optically thin environment, e.g. a corona above the disc. The 
disc structure here must be different to that assumed by disc equa- 
tions where all the energy is dissipated in the optically thick ma- 
terial. Secondly, these coronal hard X-rays can illuminate the disc, 
again changing its structure (e.g. Nayakshin, Kazanas & Kallman 
2000a). Thirdly, even without the coronal emission, the disc spec- 
trum at a given radius is not a simple blackbody as the absorption 
opacity is small at the high temperatures expected for X-ray binary 
discs (SS73). The amount of deviation from blackbody (termed a 



colour temperature correction) depends on details of the disc ver- 
tical structure, but the net result is to increase the observed peak 
'temperature' of the disc (SS73) 

A pessimistic response to this complexity is simply to aban- 
don any information the disc spectra might contain as being too dif- 
ficult to disentangle. However, some black holes can show spectra 
which are dominated by the disc component, with very little hard 
X-ray emission. These offer a way of circumventing the first two 
problems described above, so these spectra can equally be viewed 
as giving observational constraints on the disc structure, black hole 
mass and spin (see e.g. Ebisawa et al. 1991; 1994; Kubota, Mak- 
ishima & Ebisawa 2001; Kubota & Makishima 2003) 

Further motivation for this more optimistic approach is that 
the disc stability is also a sensitive indicator of how the gravita- 
tional energy is dissipated within the disc. Standard discs (a-discs) 
in which the heating is proportional to the total (gas plus radiation) 
pressure, are viscously (Lightman & Eardley 1974) and thermally 
(Shakura & Sunyaev 1976) unstable when radiation pressure dom- 
inates. In this regime the disc undergoes limit-cycle variability on 
time-scales of a few hundred seconds between a low accretion rate, 
gas pressure dominated state and a high accretion rate, advectively 
cooled slim disc (Honma, Matsumoto & Kato 1991; Szuszkiewicz 
& Miller 1997, 1998; Zampieri, Turolla & Szuszkiewicz 2001). 

Here we compile all the available Rossi X-ray Timing Explorer 
RXTE observations of several black hole binaries which show disc- 
dominated spectra to see what information can be reliably extracted 
from the data. The observed disc spectra are surprisingly simple. 
They can be well modelled by a multi-colour disc blackbody, and 
individual objects which span a large range in disc luminosity, 
idisc, generally have an observed disc temperature oc L^ c . We 
use the sample as a whole to set constraints on the colour temper- 
ature correction and limits on how this changes as a function of 
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luminosity, and show that such spectra can be used to derive black 
hole spin if the binary system parameters are well known. We use 
the data on the colour temperature correction, together with the ob- 
served stability of the light curves, to constrain the disc structure. 
We show it must be somewhat different to that predicted by the 
standard a-disc viscous heating prescription, irrespective of phe- 
nomenological modifications such as winds/jets/coronae. The data 
are also probably inconsistent with an alternative viscosity pre- 
scription in which the heating is proportional to gas pressure alone 
(Stella & Rosner 1984). 

While no current models can fully explain the observations, 
we stress that disc spectra can be a simple, useful and reliable 
probe of the accretion mechanisms. These tests will be especially 
important for the next generation of models which will calculate 
the heating ab initio from the MHD dynamo (Balbus & Hawley 
1991) rather than using the ad hoc viscosity prescriptions described 
above. 



2 DISC TEMPERATURE AND LUMINOSITY 

An optically thick, geometrically thin accretion disc around a non- 
rotating black hole has a very robust predicted spectrum. It does not 
depend on the details of the viscous heating if the gravitational en- 
ergy released at a given radius R is emitted locally as a blackbody 
of temperature T e g (R) . The total spectrum of the disc is the sum of 
these blackbodies (generally termed a multicolour disc blackbody), 
with maximum temperature T e ff, ma x occurring close to the mini- 
mum stable orbit. The luminosity of the disc is Ldi sc oc T^ s max , 
and this completely determines the disc spectrum. 

However, a given radius in the disc need not emit as a true 
blackbody. The absorption opacity at the high temperatures char- 
acteristic of X-ray binary discs is small, so scattering is important. 
In this case the disc spectrum is either a modified blackbody (with 
detailed shape depending on the vertical density and temperature 
structure of the disc) or Comptonized into a Wien peak (SS73). 
When these are taken into account, the emergent disc spectrum can 
still be approximated by a multicolour disc blackbody but of max- 
imum colour temperature, T max = / C oiTeff,m ax (Shimura & Taka- 
hara 1995; Merloni et al. 2000), where / co i is called the colour tem- 
perature correction, or spectral hardening factor, since it is always 
greater then unity. 

To find the numerical coefficient linking luminosity and tem- 
perature, we assume the disc extends down to the minimum stable 
orbit at 6R g in a pseudo-Newtonian potential (Paczynski & Wiita 
1980), 

*(*) = - GM 



Ldis 



0.583 



(i) 



where G is the gravitational constant, M is the black hole mass 
and R g = GM/c 2 is the gravitational radius. This is a simple 
and accurate approximation of the gravitational potential around 
a non-rotating black hole. The maximum disc temperature is pro- 
duced at 9.5R g with the stress-free inner boundary condi- 
tion (e.g. Novikov & Thorne 1973). Appendix A of Gierlihski et 
al. (1999) shows that the integrated disc luminosity is related to the 
maximum observed colour temperature, T max , as 



Ldis 



noG M T£ 



(2) 



Here a is the Stefan-Boltzmann constant and Co ~ 0.1067. The 
above formula can be expressed in terms of Eddington luminosity, 
i K dd = 1-48 x 10 38 (M/M Q ) erg s _1 , as 



M 



/col J V 10M © 



V 1 keV J ' 



(3) 



where k is the Boltzmann constant. 

There is considerable speculation as to whether the inner 
stress-free boundary condition is accurate. The viscosity is now 
known to be derived from an MHD dynamo, and magnetic re- 
connection can give continuous stress across the last stable orbit 
(Agol & Krolik 2000). Ironically, in this case the simple DISKBB 
model (Mitsuda et al. 1984; in the XSPEC spectral fitting package) 
is more appropriate as it does not incorporate the stress-free bound- 
ary condition (though it also has a Newtonian potential rather than 
pseudo-Newtonian). In any case, the shape of the DISKBB spec- 
trum is extremely close to that derived from the pseudo-Newtonian, 
stress-free disc model (DISKPN in XSPEC) described above, requir- 
ing only a 4 per cent temperature shift to match the spectrum at all 
energies above ~ 0.5kT max (Gierlihski et al. 1999). Thus we use 
the DISKBB model with an additional temperature correction factor 
of £ » 1.04. 

The emerging disc emission is then affected by the relativistic 
effects in strong gravitational potential of the black hole. Follow- 
ing Cunningham (1975) and Zhang, Cui & Chen (1997) we apply 
corrections g and /gr, for the observed flux and temperature, re- 
spectively: 



F' = g(i,a«) 



2-kD 2 



and 



T majc = /gr(«, a*)£T m 



(4) 



(5) 



where i is the inclination angle, and D is the distance to the source. 
The dimensionless spin is defined as a, = Jc/GM , where J is 
the angular momentum of the black hole. The primed quantities 
denote values in the observer's frame. We use the values of g and 
/gr from table 1 in Zhang et al. (1997). 



3 SAMPLE SELECTION 

We want to select disc-dominated spectra, where contribution from 
the Comptonized tail is as small as possible. We use the black hole 
sample of Done & Gierlihski (2003, hereafter DG03) as our basis. 
This consists of all the well observed RXTE black holes with low 
absorption which showed state transitions, i.e. these include soft 
state spectra. We extend the sample to include XTE J 1650-500 and 
two black hole systems which were always in the soft state when 
observed with RXTE (LMC X-l and GRS 1739-278). We also in- 
clude data from the Ginga black hole transient GS 1 124-68 (a.k.a. 
Nova Muscae), as this showed a soft (disc dominated) spectrum 
over much of its outburst (Ebisawa et al. 1994). 

To select the disc-dominated spectra from the whole sample 
we fit the data by the model consisting of the disc emission and its 
thermal Comptonization. We describe the model and fitting proce- 
dure in the next section. Ideally we would require that there was no 
Comptonized tail present. However, this sort of ultrasoft spectrum 
is not very common, so in order to have a larger sample we allow 
up to 15 per cent of the total bolometric luminosity to be present in 
the tail. This is somewhat model dependent, but the precise value 
of the coronal fraction is not important as long as it is small enough 
to have a negligible effect on the disc structure. 

Fig.Qshows the complete intrinsic colour-colour and colour- 
luminosity diagram as in DG03 for our observed sources. The lower 
curving line in the background shows the colours expected for a 
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Figure 1. The colour-colour (left) and colour-luminosity (right) diagram of the sources listed in TablefTI See DG03 for detailed description of these diagrams. 
Different symbols correspond to different fraction of the total luminosity released in the Comptonized tail above the disc spectrum: 5 per cent (black filled 
circles), 5-15 per cent (red open triangles) and > 15 per cent (cyan open circles). Only the data with disc-dominated spectra 15 per cent of Comptonized 
emission) are analysed in this paper. See electronic edition of the journal for the colour version of this figure. 



pure DISKBB spectrum with temperature increasing along the line 
from 0.5 keV (lower left) to 4.5 keV (top right); see also fig. 1 
in DG03. This line forms a lower bound on the disc-dominated 
spectra. The Comptonized fraction increases with increasing both 
colours, roughly along the diagonal. The disc-dominated data se- 
lected for further analysis are marked by black filled circles and 
red open triangles, to denote spectra where the Compton tail frac- 
tional luminosity is between 0-5 and 5-15 per cent, respectively. 
Typical X-ray spectra for the selected data are shown in panels a, 
g, h, i and j of fig. 3 in DG03. 

The limit of 15 per cent of flux in the Compton tail excludes 
all data from XTE Jl 118+480, a source well observed by RXTE 
(e.g. McClintock et al. 2001), as it has been seen only in the hard 
state. More surprisingly perhaps it also excludes Cyg X-l as al- 
though this makes a transition to the soft state, its X-ray spectra al- 
ways include a substantial fraction of Comptonized emission (see 
e.g. Gierlinski et al. 1999; Frontera et al. 2001; DG03). The inter- 
esting source GRS 1915+105 is also excluded as again the spectra 
generally show substantial Comptonized emission as well as hav- 
ing high (and probably variable) absorption which complicates the 
spectral fitting. This gives a sample of 10 black hole systems, as 
listed in TableQ In this table we also list the mass, distance, inclina- 
tion and absorption column used in the analysis, together with their 
uncertainties. There are two sets of distance estimates quoted for 
GRO J 1655-40, as the recently updated value (Mirabel et al. 2002) 
is significantly different to the previous determination. 

The intrinsic colour-colour diagram (Fig.0 shows an almost 
one-to-one match between colour and the fraction of Comptoniza- 
tion inferred from the spectrum. There is, however, one exception 
- a cluster of points from GRO J1655-40 around colours (1.25, 
0.20). They lie almost on the disc blackbody line and their spectra 
can be fitted by the disc blackbody model with temperature of ~ 
1.4 keV. However, as shown by Kubota et al. (2001), these spectra 
can be also interpreted as emission from the disc of lower temper- 
ature plus additional Comptonization in an optically thick plasma. 
The details of spectral modelling differ slightly between Kubota 



et al. (2001) and this work, in particular we use a Comptonization 
model rather than a power law to characterize the tail. However, 
the general result is the same: these particular spectra can be fit- 
ted either by a dominating disc, or by a much weaker disc plus 
strong Comptonization. In Fig. [2] we show the unabsorbed model 
components of both fits. They both give a similar \ 2 but the disc- 
dominated model (left panel in Fig.|2j yields an unphysically small 
inner disc radius (see also Kubota et al. 2001). Therefore, we select 
the strongly Comptonized model (right panel in Fig.|2j as the more 
physically realistic. Since these spectra are then dominated by the 
Comptonized emission rather than by the disc they are excluded 
from further analysis. 



4 SPECTRAL MODELLING 

For the RXTE sources we use the same data selection as DG03, i.e. 
extract one Proportional Counter Array (PCA) spectrum per one 
pointing (designated by a unique observation ID), typically giving 
exposures of a few kiloseconds. We use data from detectors 0, 2 and 
3, top layer only, except for GRO J1655-40 and XTE J1739-278, 
where detector 3 was off most of the time, so we used detectors 

and 2 only. We fit these spectra in 3-20 keV band, after adding 

1 per cent systematic error in each energy channel. For the Ginga 
observations of GS 1 124-68 we use the spectra extracted by Zycki 
et al. (1998) and fit these over the 2-20 keV bandpass with 0.5 per 
cent systematic error. 

We follow DG03 and fit the data with a model consisting 
of multicolour disc blackbody (DISKBB in XSPEC) and its ther- 
mal Comptonization (Zdziarski, Johnson & Magdziarz 1996). We 
add a Gaussian line and smeared edge in the iron Kq complex 
to approximately account for the effects of Compton reflection. 
This adequately fits all the data from all black hole (and neutron 
star) spectral states (DG03). We stress that though the model is 
highly simplified, it is appropriate for our purpose. We are not in- 
terested here in detailed spectral modelling of Comptonized emis- 
sion (which is probably non-thermal rather than thermal in the soft 
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Source Name 




D (kpc) 


% (deg) 


N H (10 22 cm -2 ) 


LMC X-3 


7(5-nr 


52(51.4-52.6)° 


60 (50-70) c 


0.07 d 


LMC A- 1 


1 A t A n <\6 

10 (4—12.0) 


JZ (Dl.4— OZ.O) 


4 J (24—04) 


n <;f 


GS 1124-68 


7 (6.4-7.6) 9 


2.8 (2.8-4)' 1 


54 (52.5-55.5)° 


0.1 6 i 


XTEJ1550-564 


10 (9.7-1 1.6) J 


5.3 (2.8-7.6) J 


72 (70.8-75.4)3 


0.65 fe 


XTEJ1650-500 


[10] 


4 (2-6) ; 


30 «40) m 


0.78" 


GRO J 1655-40 


7 (6.8-7.2)° 


3.2±0.2Por().9±0.1° 


70 (64-7 l) r 


0.8-f 


GX 339^4 


6 (2.5-10) s 


4 (2.6-5)' 


40 (20-60) s 


0.6* 


GRS 1739-278 


[10] 


8.5 (6-11)" 


[60] 


2" 


XTE J 1859+226 


[10] 


7.6 (4.6-8)" 


[60] 


0.8 1 " 


XTEJ2012+381 


[10] 


[8.5] 


[60] 


1.3* 



Table 1. The list of the sources used in this paper, together with assumed mass (M), distance (D), disc inclination (i) and absorption column (Nh ) and their 
uncertainties. The numbers in square brackets correspond to assumed values where the constraints are not known. The numbered references are as follows: [a] 
Soria et al. 2001 [b] di Benedetto 1997 [c] Cowley et al. 1983 [d] Haardt et al. 2001 [e] Hutchings et al. 1987 [f] Gierliriski, Maciolek-Niedzwiecki & Ebisawa 
2001 [g] Gelino, Harrison & McNamara 2001 [h] Shahbaz, Naylor & Charles 1997 [i] Zycki, Done & Smith 1998 [j] Orosz et al. 2002 [k] Gierliriski & Done 
2003 [1] Tomsick et al. 2003 [m] Sanchez-Fernandez et al. 2002 [n] Miller et al. 2002 [o] Shahbaz et al. 1999 [p] Hjellming & Rupen 1995 [q] Mirabel et 
al. 2002 [r] van der Hooft et al. 1998 [s] Cowley et al. 2002 [t] Zdziarski et al. 1998 [u] Greiner, Dennerl & Predehl 1996 [v] Hynes et al. 2002 [w] dal Fiume 
et al. 1999 [x] Campana et al. 2002 



state: Gierliriski et al. 1999) nor in its reflection. Since we select 
the only disc-dominated spectra, the details of the Comptoniza- 
tion/reflection modelling affect the disc results only very weakly. 
On the other hand, our Comptonization model includes a proper 
low-energy cutoff around seed photon energy, which is important. 
An extrapolation of a power law to low energies predicts much 
more flux in the tail at the seed photon energy than that given by 
a Comptonization model. This decreases the inferred disc flux, es- 
pecially for the steep tails associated with the soft states (see e.g. 
Done, Zycki & Smith 2002). 

The highest disc temperature in the sample is ~ 1.3 keV. 
The maximum power of the DISKBB model is emitted at around 
2.4fcT max . The useful observing bandpass of the PCA extends 
down to about 3 keV, below which the detector's sensitivity rapidly 
drops and systematic uncertainties become important. This means 
that the PCA can only seldom observe the peak of the disc emis- 
sion, though it can still see a large fraction of its Wien tail. For 
example, a multicolour disc with fcT max = 0.5 keV has an am- 
plitude of ~ 30 per cent of maximum at 3 keV in a vF(y) spec- 
trum. Therefore it is feasible for the PCA to effectively measure 
the observed disc temperatures down to ~ 0.5 keV. While fitting 
the spectra we limit fcT max ^ 0.4 keV (or 0.3 keV for the Ginga 
data as these have bandpass down to 2 keV) and reject observations 
for which the measured fcT max are consistent with this lower limit 
(within statistical errors) as potentially unconstrained. 

The temperature of the Comptonized component is difficult to 
constrain from the PCA data alone, particularly in the soft state, 
where the high-energy tail is weak. Leaving this parameter free 
makes deriving uncertainties of the disc parameters difficult for 
most of the soft-state spectra. Therefore, we fix it at 50 keV, which 
is well above the high energy limit of the PCA data of 20 keV. This 
makes no difference to the derived disc parameters from the disc- 
dominated spectra. 

This model fits the data very well indeed. We fit all the obser- 
vations in XSPEC, obtaining reduced % 2 jv < 1.5 (except for one 
GRO J 1655^-0 observation with strong absorption lines, see Ueda 
et al. 1998). As a result we get the disc temperature, T4 ax , and 
the unabsorbed, bolometric disc flux F ' , From this we calculate 
the luminosity in Eddington units, I/disc /£<Edd, and temperature 
corrected for relativistic effects, T max , using formulas and J5}, 




Energy (keV) Energy (keV) 

Figure 2. Unabsorbed model components fitted to the spectrum of 
GRO J1655^40 corresponding to a point in the colour-colour diagram 
(Fig.0 left panel) of colours (1.25, 0.20). The softer component (dotted 
line) represents the disc emission, the harder one (dashed line) - its Comp- 
tonization. The heavy line depicts the observed data. The two different mod- 
els presented here give a similar fit to the data. However, the disc-dominated 
model (on the left) requires an unphysically small inner disc radius. 

assuming a non-rotating black hole, a* = and mass and distance 
for each source from TableQ 



5 RESULTS 

The observed disc temperature-luminosity (hereafter T-L) rela- 
tions are presented in Figs. [3] and |4] They are roughly consistent 
with Ldisc oc T^ ax for all the sources over a broad range of lu- 
minosity. This means there can be only little change in / co i with 
accretion rate. 

In each of the figures we also plot the expected T-L relation 
from Eq. J3} for a colour-correction factor of / co i = 1.8. Obvi- 
ously, the calculated T max and -Ldisc/^Edd depend on rather un- 
certain values of the distance, mass and inclination angle (and un- 
known black hole spin). The major uncertainties from distance and 
mass are included in the figures; the cyan (or light grey in B&W) 
shadows above/below the data points correspond to the upper dis- 
tance and lower mass/lower distance and upper mass estimates, re- 
spectively. For three sources: XTE J1650-500, GRS 1739-278 and 
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Figure 3. Disc luminosity versus its maximum temperature in nine sources from Table HI (XTE J2012+361 is shown separately in Fig.|4j- The black filled 
circles and red open triangles have the same meaning as in Fig.Q The cyan (or light grey in B&W) points represent lower and upper limits on the data resulting 
from uncertainties in the distance and black hole mass (see TableHl. The lower and upper data points of GRO J 1 655^4-0 correspond to distances of 0.9 ±0.1 
(Mirabel et al. 2002) and 3.2 ± 0.2 kpc (Hjellming & Rupen 1995), respectively. The diagonal line plotted in each panel represents a relation Ldisc ?^ ax 
for a non-rotating black hole (Eq.0, calculated for the best-estimated mass (TableHl and / co i = 1.8. See electronic edition of the journal for a colour version 
of this figure. 



XTE J 1859+226, where the mass is unknown, we have assumed 
a 'canonical' black hole mass of 10 Mq with error range of 5- 
12 Mq following the mass distribution of known sources (Bailyn 
et al. 1998). In case of XTE J2012+381 there are no mass or dis- 
tance estimates whatsoever, so we plot its results in a separate figure 
(Fig.|U for representative values of M = 10 Mq and D = 10 kpc. 

We stress that though both T max and Ldisc/^Edd depend on 
uncertain parameters, for a given source its distance, mass, spin 
and inclination are constant. Therefore, the shape of the pattern 
it traces in the T-L diagram is robust. In particular, most of the 
observed sources follow a straight line in the T-L diagram (which 



is a log-log plot) corresponding to a Ldi sc <x T^ax relation. This 
can be particularly clearly seen for LMC X-3 (see also Kubota et 
al. 2001), XTE J1859+226 and XTE J2012+381. 

XTE J 1550-564 is particularly interesting, since it spans al- 
most two orders of magnitude in luminosity. It generally follows 
Ldisc ot Tn lax relation, though there is a small departure from this: 
with increasing temperature it seems to be slightly underluminous. 
This is particularly pronounced above fcT max ~ 0.9 keV where a 
break in the T-L relation can be seen. This has been also reported 
by Kubota & Makishima (2003), who noticed that the data above 
the break is consistent with Ldisc oc T^ax- They refer to this as 
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the apparently standard regime, and suggest that it is associated 
with a transition to a slim disc (Abramowicz et al. 1988). However, 
this should only occur at/above LEdd (Abramowicz et al. 1988; 
Shimura & Manmoto 2003), so it seems more likely to represent a 
subtle change in colour temperature correction. 

Other subtle deviations from a constant colour temperature 
correction can also be seen. LMC X- 1 does not show much variabil- 
ity, but over the observed small range it has an apparent anticorrela- 
tion of temperature and luminosity. XTE J1650-500 shows similar 
behaviour at the high-temperature end of the T-L diagram, where 
data points depart from the T^ax line. In both cases the observa- 
tions in question show larger fraction of Comptonized luminosity. 
They might belong to the anomalous regime (Kubota et al. 2001; 
Kubota & Makishima 2003), where Comptonization of the disc be- 
comes important. We note that the points on the low-temperature 
end of XTE J1650-500 have significant error bars and are consis- 
tent with I/di S c oc T,^ ax relation. 

The overall consistency of the disc-dominated data with 
Ldisc oc Tmax relation means there is little change in / co i with 
luminosity (see Eq. The slight departures from this relation, 
seen in GRO J1655-40 and XTE J1550-564 correspond to f col 
increasing with luminosity. This is an important observational con- 
straint on disc models. Shimura & Takahara (1995) and Merloni 
et al. (2000) both calculate the colour temperature correction ex- 
pected from a standard SS73 disc around a black hole, including 
radiative transfer through the vertical structure of the disc. Their 
predictions are shown in Fig. [5] together with the XTE J1550-564 
data. The observations are clearly inconsistent with the Merloni et 
al. (2000) results, which predict a significant decrease of / co i with 
luminosity. This shows that the vertical structure is not as assumed 
in these calculations. On the other hand, these data are consistent 
with predictions of Shimura & Takahara (1995), which are marked 
by stars in Fig. [5] 

The absolute value of / co i (as opposed to its rate of change 
with I/disc/^Edd) is somewhat more complex to constrain. At a 
given Ldisc /^Edd the temperature can be higher if the black hole is 
spinning, or if there is continuous stress at the last stable orbit as op- 
posed to the stress-free boundary condition. However, none of these 
effects will push the temperature down. For the best constrained 
source, GS 1 124-68, the data show that the maximum value of the 
colour temperature correction is ~ 2.3. For other sources it can be 
much higher. There is no stringent lower limit on / co i from the data, 
as the same T-L relation can be characterized be a black hole with 
higher spin and smaller / co i. Therefore, a conservative minimum 
value of / co i is unity, i.e. we see a true blackbody spectrum. How- 
ever, this seems entirely unlikely given the high temperature and 
associated low absorption opacity. The lowest / co i suggested in this 
case is ~ 1.6 (Shimura & Takahara 1995). Considering GS 1124— 
68 to be a representative source, we can roughly limit / co i to be 
between 1.6 and 2.3. 

Thus the use of / co i = 1.8 seems to be well justified. Many 
of the sources in Fig.|3|lie close to the T-L line for a Schwarzchild 
black hole with stress-free inner boundary condition. Again, the 
well constrained objects give the most information. A colour tem- 
perature correction of 1.8 implies that LMC X-3 and GS 1124-68 
are neither spinning nor have large amounts of energy released via 
continuous stress across the last stable orbit. Zero spin is clearly 
not a problem, but the lack of stress at the inner boundary is in con- 
flict with the numerical simulations of (non-radiative) disc structure 
(Agol & Krolik 2000). 

Several of the sources lie below the expected T-L line, most 
noticeably GRO J1655-40 (especially when D = 0.9 kpc is as- 




Figure 4. Disc luminosity versus its maximum temperature for 
XTE J20 12+361. See Fig. |5] for description of symbols and the diagonal 
line. There are no constraints for mass or distance for this source and we 
assume M = 1OM , D = 10 kpc and i = 60°. 

sumed) and GRS 1739-278. Since we have already constrained 
the colour temperature correction and stress at the inner bound- 
ary then the only parameter left to change is spin. A spinning black 
hole might also have a different colour temperature correction or 
boundary stress condition, but we assume these stay constant to 
demonstrate the effect of the black hole spin. Since the position 
of the T-L line predicted by Eq. J5} depends on the mass of the 
source, this creates another source of uncertainty. Therefore, we 
look at two sources with relatively precise mass estimates, namely 
XTE J1550-564 and GRO J1655^t0. In Fig. ©we plot their T-L 
diagrams once again, but this time assuming a maximally rotating 
black hole (a» = 0.998, with the T-L relation taken from Page 
& Thorne 1974). Clearly, XTE J1550-564 is not consistent with a 
maximal spin, while GRO J 1655-40 can be maximally rotating if 
its distance is given by the most recent determination of 0.9 kpc 
(Mirabel et al. 2002). 

All the accretion discs observed here are apparently stable 
and show no sign of the limit-cycle variability expected at high 
accretion rates. To confirm this, we examine the light curves (16- 
s resolution) corresponding to analysed spectra. They are typi- 
cally very smooth, with rms variability not exceeding 4 per cent. 
Only in the highest coronal fraction observations of LMC X-l and 
GRO J1655-40 can the rms can reach 8-10 per cent. 



6 THEORETICAL COLOUR TEMPERATURE 
CORRECTIONS 

Fig. |3| shows that in general the disc-dominated spectra have 
Ldisc oc T^ 4 . This implies a constant colour temperature correc- 
tion over a wide range of luminosity. This is somewhat unexpected 
from theoretical models of disc structure, as these predict that the 
emitted spectrum is formed under rather different conditions for 
different accretion rates, mdi sc (= ^ALjiscC 2 /LEdd, where ac- 
crestion efficiency r\ — 1/16). The first regime is at low rhdisc 
where the opacity is dominated by absorption processes and scat- 
tering is unimportant. The spectrum from a single radius is then 
very close to a simple blackbody of temperature T e s, so / co i = 1. 
What we see should be close to the disc blackbody predictions. 
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Figure 5. Disc luminosity versus its maximum temperature for XTE J 1550- 
564 taken from Fig. [3] For clarity, we don't show distance/mass uncer- 
tainties, i.e. the data points here correspond to D = 5.3 kpc and M = 10 
M© . The diagonal lines are calculated for different values of the colour- 
correction factor, / 0o ] = 1.8, 2.0, 2.2, 2.4 and 2.6 (from top to bottom). 
The overplotted thick line is the expected T-L relation from Merloni et 
al. (2000), which predicts significant decrease in / co i, with increasing ac- 
cretion rate. Contrary to this prediction, the data are consistent with f co \ re- 
maining constant (or perhaps slightly increasing) over a very wide range of 
disc luminosity. The three stars represent calculations of Shimura & Taka- 
hara (1995), which are in agreement with the data. 



At higher rhdisc the temperature rises and the absorption opac- 
ity drops, so scattering can become important. The spectrum is dis- 
torted from a true blackbody to a modified blackbody (e.g. Rybicki 
& Lightman 1979). This has a lower emissivity than a blackbody, 
but must still dissipate the same amount of gravitational energy, so 
its 'temperature' needs to be higher. Since the modified disc black- 
body spectrum depends on the absorption opacity, it is also depen- 



dent on the density structure of the plasma. The resultant spectrum 
differs depending on whether the density structure is constant with 
height (radiation pressure dominated) or Gaussian (gas pressure 
dominated: see equations 3.6 and 3.7 and fig. 2 of SS73). However, 
there is also a further effect which is that radiation can escape from 
much deeper into the disc, where the material is hotter. Both these 
effects give an observed temperature which is higher than T e g, so 

/col > 1. 

At even higher rhdi sc , the absorption opacity is so small that 
scattering dominates. Emission from deep in the disc (very hot ma- 
terial) can escape to the surface without being absorbed, but it is 
downscattered by its many interactions with electrons in the cooler 
outer layers of the disc. The spectrum is dominated by a Wien peak, 
again giving / co i > 1 (see fig. 2 of SS73). 

Merloni et al. (2000) and Shimura & Takahara (1995) show 
detailed calculations of the spectra from standard Q-viscosity ac- 
cretion discs. Both find that Compton down-scattering and the 
vertical temperature structure give rise to / cc j ~ 1.8 at high 
rhdisc 3> 0.1, where the disc is dominated by radiation pressure. 
However, their calculations differ at low rhdisc < 0.1 where the 
disc becomes gas pressure dominated (see Fig. Both Merloni 
et al. (2000) and Shimura & Takahara (1995) agree that Compton 
down-scattering is not important in these spectra, but disagree on 
the resulting colour temperature correction, giving f co \ = 2.7 and 
1.8, respectively. 

There are potential problems with both sets of calculations. 
Shimura & Takahara (1995) do not include heavy element opac- 
ity, so underestimate the true absorption in the disc, especially at 
low rhdisc. Merloni et al. (2000) do include these elements, but 
they assume that the vertical density profile is always given by that 
of a radiation pressure dominated disc i.e. is constant with height 
(SS73; Ross & Fabian 1996). However, their lowest rhdisc cal- 
culations, where the colour temperature correction becomes very 
large, give a gas pressure dominated disc. Given that these calcula- 
tions are in the modified disc blackbody regime, this would over- 
predict the temperature as the constant density discs give a larger 
colour temperature correction than the gas pressure dominated ones 
(fig. 2 of SS73). The extent of this overestimate can be estimated 
for an isothermal atmosphere from equations 3.6 and 3.7 of SS73 
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as ~ 1.66/1.2 = 1.38, which reduces the Merloni et al. (2000) 
prediction to f co i ~ 1.9. 

Thus with self-consistent vertical density structure, it appears 
that both Merloni et al. (2000) and Shimura & Takahara (1995) 
predict that the spectral hardening factor for Shakura-Sunyaev discs 
is more or less constant over the range of rhdisc explored here, as 
observed. 



7 DISC STABILITY 

Standard a-discs are viscously and thermally unstable when radia- 
tion pressure, P ra d, dominates over gas pressure, P ga s- This is be- 
cause the a-viscosity prescription assumes that the viscous heating 
rate is proportional to the total pressure P to t = Prad + Pgas- A 
small increase in mass accretion rate produces a small rise in tem- 
perature, which leads to a much larger increase in heating in the 
radiation pressure dominated regime (as P ra( j oc T 4 ) than in the 
gas pressure dominated case (where P gas oc T). Radiative cooling 
cannot keep pace with the intense heating in the radiation pressure 
dominated disc, so these a-discs are unstable at this point. The run- 
away heating can be halted when the a-disc equations are extended 
to include radial energy transport (optically thick advection). The 
advective cooling can balance even the intense P ra d heating, so 
these slim discs are stable (Abramowicz et al. 1988). 

However, while the rhdisc at which the disc becomes radiation 
pressure dominated is well defined, this merely marks the point at 
which the disc becomes locally unstable at a given radius. To trans- 
late this into global limit cycle of the inner disc requires that the 
instability can propagate to neighbouring radii, and this depends 
on details of the disc structure equations. Global numerical simula- 
tions with simple vertical disc structure show that this leads to limit 
cycle behaviour at rhdisc = 0.06 (Honma et al. 1991 ; Szuszkiewicz 
& Miller 1997; 1998; Zampieri et al. 2001), while a more sophisti- 
cated vertical structure treatment shows that the disc is globally (but 
not locally) stable up to rhdisc = 0.28 (Janiuk, Czemy & Siemigi- 
nowska 2002). 

The objects observed here span the range 0.01 < L/L-^dd 
0.5 and all of them are stable. The unusual microquasar 
GRS 1915+105 (not included in our sample) is the only source 
which shows variability which can be interpreted as this limit-cycle 
(Belloni et al. 1997). One possible explanation for this unique be- 
haviour might be that GRS 1915+105 goes to higher L/L-^dd than 
the sources considered here, i.e. that the Shakura-Sunyaev disc 
models are indeed correct in predicting that there should be a radi- 
ation pressure instability, but that they underestimate the accretion 
rate required for the onset of the limit-cycle. 

Plainly, the disc structure is different to that predicted by the 
standard a-prescription. Either there are additional cooling mech- 
anisms which act to stabilize the intense viscous heating predicted 
by these models, and/or the Q-prescription for viscous heating is 
incorrect. 

7.1 Additional Cooling Mechanisms 

The disc is cooler if some fraction, /, of the gravitational po- 
tential energy released by viscous heating is dissipated outside 
of the optically thick disc. This could also be associated with a 
mass accretion rate, rhi OS s, so that the total mass accretion rate is 
rh = rhdisc + rhi OS s. 

In the limit of rhi oss = the disc structure can be derived 
under the standard a-viscosity prescription (Svensson & Zdziarski 



1994). The disc is less luminous and cooler for a given rhdisc, or 
alternatively, for a given disc luminosity and temperature the disc 
is denser, delaying the onset of the radiation pressure instability 
(Svensson & Zdziarski 1994). The obvious candidate for this en- 
ergy loss channel is in powering a corona above the disc, since 
hard X-ray emission is generally seen from these systems. How- 
ever, we have chosen to use only spectra which have / < 0.15, 
where Svensson & Zdziarski (1994) show that the stability of the 
disc is very little changed from that of the straightforward SS73 
disc. These models do not include irradiation of the disc by this 
corona, but again the limit of 15 per cent of the power dissipated in 
the hard X-ray emission means that it has very small effect on the 
disc structure (e.g. Nayakshin et al. 2000a). 

Another candidate energy loss mechanism is to power an out- 
flowing wind or jet. Jets are commonly seen from black hole bi- 
naries except when their spectra are disc dominated (e.g. Fender 
& Kullkers 2001), ruling this out as a major energy loss route. 
Winds could be present, but line driving mechanisms which can 
power most of the mass loss inferred from the UV emitting discs 
around supermassive black holes and white dwarfs are not efficient 
in the much higher temperature X-ray binary discs (Proga & Kail- 
man 2002). A magnetically driven wind is the only remaining pos- 
sibility, in which case the results depend on the form of the wind 
(e.g. Nayakshin, Rappaport & Melia 2000b). However, such a wind 
would have to be unfeasibly powerful in order to stabilise the disc. 
Svensson & Zdzairski show that more than 95 per cent of the en- 
ergy must be dissipated in an invisible energy loss channel in order 
to produce a stable disc at radiating at Ldisc /i-Edd =0.5. 

Models in which the energy dissipation is also associated 
with some fraction of the mass accretion are no better. In the 
limit where the / = rhi OS s /rhdisc, i.e. the energy released in the 
disc and alternative loss channel are proportional to their mass ac- 
cretion rates, the disc is identical to a standard Shakura-Sunyaev 
disc with accretion rate rhdisc i.e. is unstable over the observed 
0.01 < L/LEdd < 0.5, and there is no obvious identification of 
this proposed energy/mass loss channel. A hard X-ray corona is 
still ruled out from lack of hard X-ray emission, jets are still ruled 
out from lack of radio emission, and (unobserved) winds would 
still need to be very powerful in order to stabilise the disc. How- 
ever, these models also include the energy exchange (via conduc- 
tion and radiation) between the two accreting phases in the case 
of a corona. These 'accreting corona' models calculate the verti- 
cal structure associated with an standard a-disc, and show that a 
hot skin can develop over the main body of the disc. However, this 
corona is always optically thin, and becomes negligible at the high 
L/L^dd considered here (Witt, Czerny & Zycki 1997; Rozanska & 
Czerny 2000), so has no effect on the disc stability. 

7.2 Alternative viscosity prescriptions 

The difficulties of avoiding the radiation pressure instability with 
standard a-discs make alternative prescriptions attractive. One old 
idea is to replace the a-viscosity (in which the heating is propor- 
tional to the total pressure, P to t = Pgas + Prad) with one which 
simply depends only on P ga s- Such discs (hereafter termed /3-discs) 
are stable even when dominated by radiation pressure (Lightman & 
Eardley 1974; Sakimoto & Coroniti 1981; Stella & Rosner 1984), 
since the heating is always oc T. 

While such models trivially solve the disc stability, they give 
a very different disc structure and hence predict different behaviour 
for colour temperature correction as a function of accretion rate. 
For low rhdisc, where gas pressure dominates, then the (3- and a- 
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disc structures are identical, hence should have / co i ~ 1.8 (see 
Section[6j. However, at high rhdisc, where radiation pressure dom- 
inates, the density of the /3-disc is much larger than that of the stan- 
dard a-disc (Stella & Rosner 1984), though the disc thickness is un- 
changed (Stella & Rosner 1984). Compton down-scattering is not 
important in these spectra at any rhdi sc (Nannurelli & Stella 1989), 
so the disc spectrum is in the modified blackbody regime. Its shape 
at any radius depends on how the density and temperature change 
as a function of height in the disc (SS73). While full calculations 
have yet to be done, we can use the a-disc calculations in a compa- 
rable regime (radiation pressure dominated, but where the emission 
is not strongly Comptonized) to estimate the resulting colour tem- 
perature corrections. This corresponds to the lowest rhdisc of the 
Merloni et al. (2000) calculations, i.e. to / co i ~ 2.7. Thus it seems 
that the /3-disc viscosity predicts / oc j increasing from ~ 1.8 at low 
rhdisc to ~ 2.7 at high rhdisc- This is strongly at odds with what 
is observed, leading us to conclude that the aP gas heating law is 
ruled out observationally (see also Janiuk et al. 2001). 

Thus none of the ad hoc viscosity prescriptions appears able 
to give a disc structure which matches both the observed colour 
temperature corrections and stability. The answer probably lies in 
the next generation of accretion disc models which will superseded 
these ad hoc viscosity prescriptions. The physical mechanism for 
angular momentum transport is now known to be the magneto- 
rotational instability (MRI: Balbus & Hawley 1991). Numerical 
MHD simulations of the global disc structure of gas pressure dom- 
inated discs show that the a-prescription applies approximately in 
the disc body (with vertically integrated (a) = 0.1). However, 
these also show strong gradients in the ratio of viscous heating to 
pressure in both radial and vertical directions, and large fluctuations 
in the heating as a function of time at any point (e.g. Hawley & Bal- 
bus 2002). For the bright discs considered here, radiation pressure 
strongly affects the dynamics, and the mechanics of the MRI in- 
duced turbulence changes (Agol & Krolik 1998; Blaes & Socrates 
2001). The coupled radiation and magneto-hydrodynamics is a dif- 
ficult and very computationally intensive problem, but preliminary 
results indicate that the average heating is unlikely to be propor- 
tional to total pressure, as in the standard a-disc (Turner, Stone & 
Sano 2002; Turner et al. 2003). Even if it were, the MRI induced 
turbulence and associated convection may lead to enhanced cooling 
rates (Turner et al. 2002; 2003) 

Self-consistent disc models open up the possibility of being 
able to calculate the viscous heating and convective turbulent cool- 
ing associated with the MRI. When these models are sufficiently 
robust to predict the vertical disc structure and associated radiation 
spectra, then the data presented here offer a crucial test. 



8 BLACK HOLE SPIN 

'Black holes have no hair' : they should be purely characterised by 
their mass and spin (any charge would always quickly be neutral- 
ized by accretion). The stellar remnant black holes all have very 
similar masses, so the only free parameter to give observable dif- 
ferences at the same (steady state) accretion rate is spin. Prograde 
spin drags the last stable orbit closer to the black hole, resulting in 
higher luminosity and temperature for a given accretion rate onto a 
spinning black hole compared to a Schwarzschild one. This shifts 
the expected T-L relation to the right (Fig.|5J, giving an observ- 
able determination of the black hole spin (e.g. Zhang et al. 1997; 
Gierlinski et al. 2001). Retrograde spin has the opposite effect, but 
is highly unlikely since the binary stars formed together from a 



single cloud, so spin and orbital angular momentum should be cou- 
pled. 

Plainly, most of the black holes in this sample are consistent 
with Schwarzchild, lying on or close to the zero-spin T-L relation. 
LMC X-l, LMC X-3 and GS 1 124-68, where the system parame- 
ters are tightly constrained, have very little room for even moderate 
spin. Conversely GRO J1655-40 and GRS 1739-278 have signifi- 
cant spin. 

All these (apart from GRS 1739-278) were known previously: 
GS 1124-68 (Cui et al. 1997), GRO J1655— 40 (Cui et al. 1997; 
Gierlinski et al. 2001; Kubota et al. 2001), LMC X-3 (Cui et 
al. 1997; Kubota et al. 2001), LMC X-l (Gierlinski et al. 2001). 
What is new here is the size of the sample, which shows clearly 
that a constant colour temperature correction is a fairly good ap- 
proximation. The rest of the objects have rather large uncertainties, 
and are compatible with moderate (but not extreme) spin. How- 
ever, the fact that the best estimates for the system parameters put 
XTEJ1550-564, GX339-4, XTEJ1650-500 andXTE J1859+228 
so close to the Schwarzchild line is certainly suggestive, if not com- 
pelling. 

Better system parameters from optical/infrared studies of the 
systems will reduce the uncertainties, and provide a clear deter- 
mination of the black hole spin. This is especially important (but 
for different reasons) for two of the systems here. XTE J 1550-564 
has a superluminal radio jet (Hannikainen et al. 2001; Corbel et 
al. 2002), so confirmation of its Schwarzchild nature would provide 
a conclusive counterexample to the paradigm in which relativistic 
jets are triggered by black hole spin (e.g. Moderski, Sikora & La- 
sota 1998). XTE J 1650-500 is claimed to have maximal spin from 
detailed fitting of the iron line profile in the X-ray spectrum (Miller 
et al. 2002). Currently, with no good mass estimate our data can- 
not constrain its spin, though in future it might give an independent 
measurement. Here the issues are even more emotive, as the line 
profile in XTE J 1650-500 is similar to that of the AGN MCG-6- 
30-15 (Wilms et al. 2001), the defining example of a spinning black 
hole. We stress that the ultrasoft X-ray binary spectra offer an inde- 
pendent constraint on the spin, and hence give a check on iron line 
modelling. 



9 CONCLUSIONS 

We observe discs emitting in the luminosity range 0.01 < 
^disc/iEdd ^ 0.5. Individual objects which span a large range 
in I/disc/^Edd are generally consistent with Ldisc oc T 4 , imply- 
ing a constant colour temperature correction. In LMC X-3 and 
GS 1124-68, the distance, mass and inclination of the black hole 
are well constrained, giving / co j ~ 1.8, and, by implication, re- 
quire a Schwarzchild black hole. Conversely, for GRO J 1655-40 
and GRS 1739-278 the disc temperature requires moderate-to-high 
spin, i.e. a Kerr geometry. For black hole binaries which show disc- 
dominated spectra this gives an alternative way to check the inner 
disc geometry derived from modelling the relativistic smearing of 
the iron line and reflected spectra. 

The observed range of 0.01 < I/disc/iEdd < 0.5 is where 
standard Shakura-Sunyaev models predict that the inner disc should 
be unstable due to the dominance of radiation pressure. The mod- 
els predict limit-cycle behaviour in this range, yet observations 
show that the light curves are characterized by very little variabil- 
ity. These observations strongly require that the disc is stable in this 
range of L disc /L E dd- 

Stability could be retrieved in the Shakura-Sunyaev disc mod- 
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els if the disc remains gas pressure dominated. This can be arranged 
if there are other cooling mechanisms (corona, jet, winds) which 
give a much denser disc for a given mdi sc - However, since we 
choose only disc-dominated spectra, the coronal emission is neg- 
ligible. These spectra are also always observed to have strongly 
quenched radio emission, so limiting the jet emission, while winds 
would have to be unfeasibly powerful. Hence there is only very 
limited potential for avoiding the conclusion that radiation pressure 
does dominate in our high I/disc/^Edd discs. 

The only other known way to circumvent the radiation pres- 
sure instability is to change the viscosity prescription so that the 
disc heating is proportional to P gas rather than the standard a- 
disc prescription where it is proportional to P to t = Pgas + Pad- 
This avoids the tremendous increase in heating rate so the disc can 
remain stable even when dominated by radiation pressure. While 
proper calculations have yet to be made, it appears that such discs 
should have a colour temperature correction which changes ~1.8 
to ~2.7 as mdisc increases. Conversely, the standard a-discs pre- 
dict that the colour temperature correction should remain stable at 
~1.8 for the range of rhdisc considered here. The observation that 
the colour temperature correction remains stable over a wide range 
in Ldisc/ ^Edd appears to rule out the alternative viscosity prescrip- 
tion, and strongly supports the vertical structure as predicted by a 
standard Q-disc. 

This leads to an impasse. Observations show that the disc 
is most likely radiation pressure dominated, with standard a- 
viscosity. Such discs should undergo limit cycle behaviour, yet are 
observed to be stable. There are strong indications that the solution 
to this lies in the self-consistent heating and cooling generated by 
the MRI turbulence which is the physical source of the viscosity. 
While coupled radiative magneto-hydrodynamic calculations can- 
not yet predict the disc vertical structure in a robust way, this will 
probably become possible in the next few years. 

Irrespective of the outcome of detailed calculations, it is clear 
that accretion disc spectra can constrain disc models. While spectra 
which contain a substantial Compton component at higher energies 
are indeed ambiguous, we challenge the Merloni et al. (2000) asser- 
tion that fitting the accretion disc spectrum gives results which are 
not in general directly related to the actual disc parameters. As also 
shown by Ebisawa et al. (1994), Kubota et al. (2001) and Kubota 
& Makishima (2003), the disc dominated spectra can be a reliable 
guide to the accretion disc structure. 
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